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We investigate 't Hooft-Mandelstam monopoles in QCD in the presence of a single classical 
instanton configuration. The solution to the Maximal Abelian projection is found to be a circu- 
lar monopole trajectory with radius R centered on the instanton. At zero loop radius, there is a 
marginally stable (or flat) direction for loop formation to 0(R 4 logR). We argue that loops will form, 
in the semi-classical limit, due to small perturbations such as the dipole interaction between instan- 
ton anti-instanton pairs. As the instanton gas becomes a liquid, the percolation of the monopole 
■ loops may therefore provide a semi-classical precursor to the confinement mechanism. 
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I. INTRODUCTION 



Both instantons and magnetic monopoles are thought to play an important role in the Euclidean description of the 
Yang-Mills vacuum. In QCD instantons provide the solution to the U(l) problem Q and a vacuum composed of a 
liquid of instantons ^] can explain chiral symmetry breaking and the low mass spectrum, at least at a qualitative 
level HQ- Nonetheless, it is generally believed, as conjectured by Mandelstam || and 't Hooft that magnetic 
monopoles, not instantons are essential for confinement in QCD and related gauge theories. For example, compact 
QED with a lattice cutoff is known to be exactly dual to a Coulomb gas of monopoles, which upon condensation causes 
confinement via a dual Meisner effect 0. Similarly there is evidence for the role of monopole condensation in the 
3-d Yang-Mills-Higgs (or Georgi-Glashow) model ||] and more recently in 4-d N = 2 supersymmetric Yang-Mills 
! theory ||. 

Consequently, we are faced with a dilemma of two competing pictures of the QCD vacuum (or more accurately the 
Euclidean equilibrium phase) as either a coherent ensemble of instantons or a condensate of magnetic monopoles. In 
this paper, we show that these two pictures may in fact be two descriptions for the same phase. To establish a precise 
| and indisputable link between the monopole trajectories and instantons, we define the monopole current in a field 
configuration via 't Hooft's Maximal Abelian (MA) projection and demonstrate that in the background of a single 
instanton the MA projection leads to a circular monopole current loop of radius R centered on an isolated instanton 
of "width" p. We also begin the study of monopole trajectories in the presence of interacting instanton pairs. 

Let us review briefly the attempts to identify monopole currents in QCD. (The reader is referred to Polikarpov 
and references therein for more details.) It is well known that the only topological stable solutions to Euclidean Yang- 
Mills theory are the multi-instanton solutions with topological charge Q. On the other hand, the 't Hooft-Polyakov 
monopole solution to 3-d Yang-Mills-Higgs theory can also be viewed as a static solution to the pure Yang-Mills 
Euclidean field equations in 4-d, where Aq plays exactly the same role as the adjoint Higgs fieldQ in the BPS limit. 
Of course, since the classical theory is now scale invariant the monopole mass M is an arbitrary parameter. Quantum 
effects will be necessary to set the scale. Subsequently Rossi has pointed out that this monopole configuration 
can be identified with (i.e., has the same field configuration as) an infinite sequence of instantons equally spaced by 
M~ 1 8w 2 /e 2 on the time axis in the limit p — ► oo. 

Turning now to the MA projection, Chernodub and Gubarev |l4[ have recently shown that a straight line monopole 
trajectory through the center of an instanton (or pair of instantons) satisfied the MA gauge condition. However if 
one looks at the MA gauge fixing functional, G — j J d 4 xA+(x)A~(x), A^ = A 1 ^ ± iA^, for the Chernodub-Gubarev 
solution, one notices that it diverges at large distances. Since the proper definition of the MA projection is not just 



1 Essentially the same use is made of the adjoint field Ao in the MA projection that diagonalizes the Polyakov loop as discussed 
by Suganuma et al. jl2| , but we do not pursue this approach further since it does not provide a Lorentz invariant definition of 
monopole variables. 
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a stationary point but a minimum of this functional, this solution is still not a very satisfactory linkage between an 
instanton and the monopole. The basic problem is that the instanton is restricted to a space-time region whereas 
the static monopole solution has an infinite trajectory. In contrast we have found that our monopole loop solution 
to the MA projection, which is localized at the instanton, does give a finite value to G, which drops gradually to an 
absolute minimum as the radius decreases: G ~ An 2 p 2 [\ + 4.6(R/ p) log(p/ R)]. We then argue that this solution is 
easily stabilized by small Gaussian perturbations or nearby interactions with anti-instantons. We feel this is the first 
rather precise connection between the instanton and monopole trajectories. 

In addition to this continuum analysis, there has been some recent numerical evidence in lattice simulations pre- 
sented by Hart and Teper |[{|, by Bornyakov and Schierholz Jlq ] , and by Markum et al. and Thurner et al. fTij] 
that shows that monopole loops are correlated with multi-instantons configurations. These numerical results also 
indicate that a correlation between monopole condensation and the formation of a dense liquid of instantons are a 
property of the quantum vacuum. However our immediate goal is a more modest one: establishing a clear kinematical 
connection between monopoles and instantons in the semi-classical limit. While this limited exercise is unable to 
show that the dynamics of the full quantum theory is in some sense approximated by a monopole-instanton correlated 
background and that confinement results from the percolation of these objects, we feel that clarifying the kinematics 
is an important first step in bringing these two pictures together. 

To go beyond the analytical study, we have also undertaken a numerical MA projection in the background of an 
isolated instanton by minimizing the functional G discretized on a 4-d hypercubic lattice. The only solutions we have 
found were identical to our analytical monopole loop solution albeit stabilized to a fixed radius (in lattice units) by 
lattice artifacts. We have also begun to study the monopole loops in the background of an interacting instanton 
anti-instanton (I-A) pair and found that the individual loops appear to be stabilized by the interaction, if the I-A 
molecule is oriented so that the dipole interaction is attractive. At a critical separation of the I-A pair, the individual 
monopole loops fuse into a single loop. Details of the multi-instanton study are left to a future publication Jl8| . 
Nonetheless on the basis of our preliminary analysis, we conjecture that the I-A interaction provides a semi-classical 
mechanism whereby the monopoles are "liberated" from the individual instantons, and which may therefore represent 
a precursor to condensation and confinement. 

The organization of the paper is as follows. We begin by noting that the Maximal Abelian projection, which is 
usually presented as a gauge fixing procedure is equivalent to the introduction of an "auxiliary Higgs" field $ = (j> ■ t, 
which is determined by minimizing the Higgs kinetic term, 

G = \J d*x[D„(A)$} 2 , (1) 

(The Higgs field is related to the usual gauge rotation by $ = This approach has several advantages in terms 

of analysis and a clearer relation to the standard discussion of monopole topology. (See Appendix A for details.) 

Next in Section 2, we solve the stationarity equations for the Higgs field (or gauge rotation in the SU(2)/U{\) coset 
space) , to find the monopole loop solution for fixed radius R. Monopole loops in Abelian gauge theories are reviewed 
to guide the construction. 

In Section 3, we consider the full manifold of solutions and their collective co-ordinates. We note that there is 
an interesting correlation between the average over the orientations of all monopole loops at fixed R and the isospin 
orientation of the instanton given by 

< Npy >= |[(^,3^i/,4 - <5 m ,3<5ia4) _ (<W<^,2 _ <W<W)] ( 2 ) 

where iV M „ is the skew symmetric unit tensor defining the plane of the loop in 4-d. 

In Section 4, we report on our numerical solutions leading to a preliminary picture of the role of instanton anti- 
instanton (I-A) interactions for loop stabilization at large separation and for the I-A pair loop fusion and percolation 
at small separations. 

Finally in Section 5, we discuss our results and suggest future lines of investigation. We comment on the naturalness 
of the instanton-monopole connection. The loop "VEV" naturally breaks the SO(A) group down to the coset direct 
product SUl(2)/Ul(1) x SUr{2)/Ur{V) in close analogy to isospin breaking in the 't Hooft-Polyakov monopole 
construction. Also we note that the Abelian projection for the instanton in the singular gauge spreads the singularity 
at the origin to the trajectory of the monopole loop, directly relating the instanton's topological charge Q to the 
monopole charge g by the identity Q — eg/4n. In summary we conclude that monopole loops are peculiarly well 
matched to the instanton, leading us to hope that there is a deeper connection in confining gauge theories transcending 
our particular construction. 
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II. MONOPOLE LOOP SOLUTION IN THE INSTANTON BACKGROUND 



The Maximal Abelian projection was introduced by 't Hooft || in order to define magnetic monopole coordinates by 
a partial gauge fixing procedure that leaves the maximal Abelian subgroup free. In analogy with the Higgs system, 't 
Hooft suggested introducing an adjoint field X, which by the gauge transformation, X — > X n = flXW, is diagonalized 
to fix the gauge in the coset space SU(N)/U(1) N ~ 1 up to U(l) factors. Exceptional space-time trajectories where 
X has degenerate eigenvalues represent the Abelian monopole configurations. Various examples for X were suggested 
such as F12 or the (untraced) Polyakov loop, etc. 

Subsequently a particular Lorentz covariant gauge has proven to best correlate the monopoles condensate with 
confinement in lattice simulations Mj. This gauge is now referred to as the Maximal Abelian (MA) gauge. For 
SU(2), it is expressed as the minimization of the functional, 

G[A,] = lJd i xA+(x)A-(x), (3) 

A^ = A* ± 2A 2 . In differential form the MA gauge condition becomes 

A ± (x) = (d fl ±ieAl)A^ = 0, (4) 

with the additional stipulation to avoid Gribov copies that one should find the solutions corresponding to the global 
minimum of G. (We will sometime refer to the stationarity condition (|J) by itself as the "differential" MA projection.) 

Before presenting the details of our monopole loop solution, it is possible (at least with hindsight) to see why 
magnetic loops might appear in an instanton background. Consider the field configuration for a single instanton p(H22] 
in the singular gauge, 

An = \A^ • t = §T Q 7^AZn(l + p 2 /x 2 ), (5) 

where the anti-self-dual 't Hooft symbol is fj® = (i/2)Tr(T a T u T^) and t m = (i,f) with Latin indices restricted to the 
vector components. First one notices that this field configuration already satisfies the differential MA (as well as the 
Lorentz) gauge condition. Also it gives a finite value, 

G[A„] = 2e 2 J d*x x2{ / +p2)2 = 4. 2 e y , (6) 

to the MA gauge fixing functional — the dominant contribution to G coming from the core of the instanton, < x < p. 
On the other hand, the instanton in the non-singular gauge also satisfies the differential Eq. ([|) for the MA projection, 
but now the functional G diverges logarithmically at large distances, 

G[A^)] = 2e 2 |^ I - J ^- ¥ -oo. (7) 

Nonetheless, the non-singular gauge does reduce the contribution to G in the core of the instanton. Consequently, 
comparison between the singular and non-singular gauges suggests that it might be possible to find an intermediate 
solution which minimizes G by a gauge transformation such that the central region, inside some radius (x < R), 
is converted to the non-singular gauge, while the large distance behavior is unchanged from the singular gauge 
configuration. In essence this will turn out to be the way we construct the monopole loop solution. 

In addition, by a general argument one can understand why this construction might lead to monopole singularities. 
No local gauge transformation on the instanton can change the total topological charge. However if we express the 
topological charge for the instanton in the singular gauge via Gauss' law, Q — J d A xd p K p , in terms of the flux for the 
gauge variant current, 

K» = ^e^pxTr[A v {d p A x - §A p A A )] , 

one knows that the entire contribution comes from the singular point at x = 0. Any gauge rotation that removes this 
singularity but does not change the gauge at infinity must replace it with another singularity at some finite distance. 
As we will explain in more detail in the conclusion, our monopole loop current is the source of this singularity and 
its magnetic charge can therefore be directly related to the topological charge of the instanton. 
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A. General Equations 



In searching for solutions to the MA projection, we have found it more convenient to use a gauge covariant formula- 
tion in terms of an auxiliary Higgs-like field, Q(x) = SVr^n, instead of working directly with the gauge transformation 
f2(x) itself (see Appendix for details). In terms of the isovector field (f>, $ = <f> ■ f , the functional G becomes 



G=hld 4 x{ 1 1 [D^A)$f + V($ 2 )}, (8) 



where D^ a = d^ a + ee a ^Af l (jf 1 and the potential is V((j> 2 ) = a(cf 2 - 1) with a La grange multiplier a. (One can see 
by inspection that Eq. (^) is identical to Eq. (||), after the substitution $ = il^T^il.) It is also natural to generalize 
the MA projection with a quartic Higgs potential, 

v^) = i\(^-vr, (9) 

with &(x) now given by 4>(x) ■ t/\4>(x)\. We will refer to this more general from as the "Higgs MA projection". The 
conventional MA projection is recovered in the limit A — > oo at fixed VEV v in which the Higgs mass scale mn —* oo 
as well. Not only does the covariant formulation have numerous technical advantages, it shows that the MA projection 
need not be seen as a gauge fixing prescription but instead as a gauge covariant method of identifying the appropriate 
magnetic variables. 

The general problem is to construct solutions to the differential form of the MA projection, 

D fl (A) 2 $-V <l> V^)=0, (10) 

in the background of a single instanton. For the standard MA projection, this yields a linear PDE for <fi(x), 

d 2 J+ 2A„ x A^ x (A M x$)-(t$=0, (11) 

subject to the constraint <fi 2 — 1. For an instanton in the singular gauge (^), the potential also can be written as 
A-ti = -■ix 2 f(x)g' i (x)d^g(x), where g*(x) = -ix^/\x\, t m = (f,i), 

™ = ?1*T?> (12) 

and we have also set the gauge coupling e = 1. 

Solving this equation is very similar to the problem, first considered by 't Hooft Jl[, for computing the Gaussian 
fluctuations of gauge, fermionic and Higgs fields in the background of an instanton. In studying self-dual solutions and 
their supersymmetric extensions, it is useful to introduce tensors that reflect the chiral SUl(2) x SUr(2) decomposition 

of the 0(4) Euclidean Lorentz group, which suggests the conformal coordinates, u = x + iy = ue li P and v = z + it = 

ve^ , that enter into the bispinor, 

, * _ / z + it x - iy \ _ fv u* \ , 

x a - x A T ») a p - y x + w _ z+it j-y u -v* ) 

for the singular gauge transformation, g* (%) a g = —ix a &/\x\. This notation, as we will see in Sec. 3, is also convenient 
for exhibiting the Lorentz symmetries of the monopole solution. 

It follows from the property of fft y that isospin components of A^ are orthogonal, i.e., A^A^ ex 8 a p. As a 

consequence, the last two terms in Eq. ( pd[ ) both point along <j> in the isospin space. Let us parameterize this "radial 
direction" at each a; by a unit- vector, 

(j)(x) = (sin (3 cos a, sin (3 sin a, cos (3) (14) 

and introduce two unit tangent-vectors orthogonal to <j>(x), a(x) = (— sin a, cos a, 0) and (3(x) = 
(cos P cos a, cos (3 sin a, — sin (3). By projecting Eq. ( |TT| ) onto the a(x) and (3(x) axes, the two independent differ- 
ential equations for the MA gauge are 

sin/3 &la + 2cos/3 (d„a){d^) = -2sin/3(0 • (15) 
8lf3 - \ sin(2/3) (d^a) 2 = 2 sin (3(0 ■ (0„a), (16) 



4 



respectively. (The more general equations for the Higgs MA projection for arbitrary A are given in Eqs. ( A15| - A17).) 
Our task is to find solutions to these equations. 

For the explicit monopole content of a solution, it is instr ucti ve to return to the active transformation, Eq. ( |Al| ), 
and to decompose the transformed field into two parts, Eq. (A2), A^(x) — Ap(x) + M^(x) where 

A^x) = n(x)A^x)rt(x) , M^x) = ^n(x)d^(x) . 
The gauge rotation in terms of three Euler's angles is 

n(x) = e luJ ^ = e HT3/2^T 2 /2 e iaT 3 /2_ (17) 

One can readily verify that this is consistent with Eq. ( [l4|) and, owing to the residual U(l) invariance, <p{x) is 
independent of the angle 7. Focusing on the Abelian field a ll (x) = Tr(r 3 A^) — A^ + M^, one observes that the first 

term A^ = cf) ■ A^ provides the source for our PDE's, Eqs. ( |l5| ) and ([l6|). The induced term, 

M*(x) = -~[cos(3(x)d„a(x)+d^(x)], (18) 

will give rise to topological current for the monopoles, Eq. (|A5|). In terms of the Euler's angles, the magnetic current 
is explicitly given as, 

k n( x ) = T~~ € ^pad v [d a a{x)dp cos (3(x)}. (19) 
Zne 

Note that is not specified uniquely due to the residual U(l) gauge symmetry. Generally, we will choose j(x) = 
— a(x) so that the Dirac sheet associated with the monopole loop solutions can be oriented conveniently and so that 

Q.(x) = cos(/3/2) + i sin(/3/2) [(cos a)n - (sin a)r 2 ] (20) 

obeys the boundary condition — > 1 as (3 — > 0. 



B. Monopole Loop in 17(1) Gauge Theory 

To construct monopole loop solutions, we have found it helpful to work "backward" from specific examples of 
U(l) monopole currents. First, recall Dirac's construction for a static monopole of magnetic charge g at the origin 
x = y = z — 0, with magnetic current, k fl (x) — gS^'fflS^. The field is given by 



JJ_ 

4-7T 



[l-costf]d^, (21) 



in spherical coordinates, < t? = tan 1 (y / x 2 + y 2 / z) < tt and <p — tan 1 (y/x). The singular behavior due to 
dfj_(p = if/ \/ x 2 + y 2 as \/ x 2 + y 2 — ► corresponds to having placed a Dirac string along the negative z-axis, which in 
4-d leads to a Dirac sheet in the left-half of the 3-4 plane. The Dirac sheet can be moved by a gauge transformation. 

This construction is an example of a general formalism for writing down a solution with a Dirac sheet attached to 
a closed Euclidean monopole current loop ^3|. Consider a unit monopole trajectory x M = y^r) and its associated 
magnetic current k^{x) = J d,Td T y lll (T)5 4 (x — y). We can attach to it a Dirac sheet described by t/ m (t, a), where 
%n = Hfj.(T,0)- The associated vector potential for the monopole loop is 

/ \ 1 f j4 i^vG v ^(x') 

a >W = (22) 

where = ^^G^and G^x) = J dTda{d T y ll d IJ y lJ - d T y v d (T y IJ }8 i {x - y(r,a)). 

We now apply Eq. ( J22| ) to a monopole of charge g moving in a closed loop of radius R in the 3-4 plane centered at 
the origin. The Dirac sheet (or solenoid world sheet) is chosen to run across the loop (see Fig[l]), parameterized by 
yi=?/2 = 0, 7/3 = R(l — a) COS7TT, 7/4 = i?(l — a) sin7TT, where —1 < r < 1, and < a < 1. It is convenient to use 
the conformal coordinates introduced earlier, u = x + iy = ue 1 ^ and v = z + it = ve^. The magnetic current, 

K{x) = -jt-8{x) S(y)S(v - R)$, (23) 
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results in the potential 



tt = T Z [l-«»(e + + 0-)]d ll <p i (24) 



4m 

where 6± = tan _1 [u/ (v ± R)]. (See Fig [j] for a geometric interpretation for 9±.) 
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FIG. 1. (a) The Dirac sheet for a monopole loop in QED. (b) 8± variables for Dirac loops solution. 

As for the case of a single monopole, d^ip = (p/u is singular as one approaches the 3-4 plane, u — > 0. For v > R, 
6± — > 0, the pre-factor vanishes, and this singularity is absent. However, for < v < R, 9 + — > 0, 6L — > 7T, the 
pre-factor approaches 2 and a Dirac sheet is present. We are now ready to generalize this construction to the Maximal 
Abelian projection in the background of an instanton. 

C. Monopole Loop in an Instanton Background 



Returning to the differential MA conditions, Eq. (|15|) and Eq. (|16|), an obvious solution corresponds to the gauge 
rotation from the singular to the non-singular gauge, f2 = = {xa + ix ■ t)/\x\. In the Euler parameterization, this 
corresponds to a = ip — ip, 7 — it — (</? + ip), and (3 = 29 = 2tan _1 (u/w). However, on surfaces where d^a is singular 
(u = or v = 0), cos (3 = so that no magnetic monopole is formed.^ 

With the gauge choice 7 = —a, the induced potential, M^, becomes 

M u(x) = -[1 - cos/3(a;)]a /1 a(ar). (25) 



Note the similarity between this expression and the monopole solutions for QED, Eq. (|2J). The Dirac sheet is formed 
on the surface where d^a becomes singular and 0{x) = ±tt. The magnetic current is on the boundary of the Dirac 
sheet where (3 has a discontinuity. 

Consider the ansatz where a and (3 are only functions of (cp, and (u, v) respectively, i.e., a(ip,ip) and (3(u,v). It 
can be shown that the first of the two MA conditions, Eq. (Ea), is solved under this ansatz by 



2 The solution J1J] given by Chernodub and Gubarev corresponds to [3 = 1? and a — tp, where ■& and (p are the polar and 
azimuthal angles for the spatial three vector x. This is the standard static "hedgehog" configuration for <f> a = 1° /\x\, which 
gives rise to the 't Hooft-Polyakov monopole with charge Aiv/e at x = 0. This solution leads to a divergence value for the gauge 
fixing functional G. 



G 



a = ip — ip. 



(26) 



It follows from the QED example that a Dirac sheet can be present either in the 1-2 plane, (v — > 0), or the 3-4 plane, 
(u — > 0), or both. It is sufficient for us to seek solutions where the monopole loops are oriented in the 3-4 plane. Other 
orientations, including that oriented in the 1-2 plane, can be obtained by performing appropriate 0(4) rotations, as 
we demonstrate in Sec. 3. 

Instead of u and v, we convert to x and 8, where x 2 = u 2 + v 2 and 6 = tan _1 (u/w), < 8 < ir/2. (See Section 4 
for further discussion of the ansatz in the u-v coordinate system.) Now Eq. ([l6]) becomes 

±-d x (x 3 d x [3) + 1 de(sm20de/3) - „ 2 + 4/(s) sin 2 /?(cot (3 - cot 26) = 0. (27) 
x A rsin2S x 2 sin (26) 

For G to remain finite, one must impose the boundary condition sin/3 = at 8 = and 8 = 7r/2. There remains the 
freedom for (3 to take on values which are integral multiples of 7r. 

Recall that we are seeking solutions with A^j(x) approaching the behavior of the non-singular gauge at the origin, 
and approaching the behavior of the singular gauge at infinity. We therefore expect a solution with (3 ~ 28 (modi:) for 
x small and with (3 — > (mod n) at infinity. A monopole loop lying in the 3-4 plane corresponds to a solution where (3 
has a discontinuity on the positive w-axis. To be precise, we seek solutions where (i) (3(x,tt/2) = 0, (ii) P(x,0) = — tt 
for < x < R, and (3(x, 0) = for R < x < oo, and (iii) (3(x, 8) — > for x — > oo. 

An analytic solution for (3(x, 8) in the limit of small monopole size, Rj p — > can be found; this will be discussed in 
Section 3. For general R/p, the solution to Eq. ( p7| ) can be obtained numerically. Using a simple relaxation method 
for G (described in Section 4) , very accurate solution for (3 can be found, as illustrated in Fig ^. Independent of the 
detailed form of the solution, with a = ip — ip, the discontinuity in f3(x, 6 = 0) at x = R leads to a magnetic current, 

= - S(x 1 )d(x 2 )S(v - (28) 
e 

where we have re-introduced the charge factor. Therefore our solution corresponds to a loop of monopole with a 
magnetic charge g — An/e in the instanton background. 
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FIG. 2. Solution to f3(u,v) for a monopole loop with R/p = 0.81 and p = 1. 

In fact we have explored the more general case of a Higgs potential V(4>) — X((j> 2 — l) 2 , where 4> is the magnitude 
of the Higgs field cj). The change in the value of G is almost independent of the parameter A. It is now clear that 
our general monopole loop solutions yield a family of solutions which all satisfy the local MA condition and which 
interpolate between the singular gauge and the non-singular gauge for an instanton. For < x << R, f2(x) ~ g (x) 
up to a U(l) rotation; the effect of the gauge transformation is to remove the singularity at x = 0. Since (3 — > 
for x — ► oo, it follows that, with 7 = —a, Q. — ► I rapidly and the large-x behavior of the gauge fields is unchanged. 
The transition between "small" x and "large" x is marked by the presence of a monopole loop of radius R in the 3-4 
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plane. As R goes to oo, the monopole loop is pushed to infinity, leaving behind a gauge field configuration which is 
the instanton in the non-singular gauge. 

There is in fact a very large manifold of solutions. In addition to the loops described above, one can clearly impose 
the boundary conditions for (3 to jump by ±7r on either the 9 = axis or the — 7r/2 axis or both, at arbitrary 
radii Ri. So long as the final value of (3 as x 2 — ► oo is set to be the same multiple of it on both boundaries, no 
"topological charge has been pushed to the spatial infinity and the functional G will still be finite. This will give a 
series of concentric loops with increasing radii and magnetic charge either g or —g. In addition as describe in Sec 4, 
the application of Lorentz invariance to these solutions will re-orient these loops in 4-d space. 

III. COLLECTIVE COORDINATES AND STABILITY OF MONOPOLE LOOP 

In this section, we first address the question concerning the degeneracy of the monopole loop solutions constrained to 
a fixed radius R in the instanton background. This analysis closely parallels the discussion of the collective coordinates 
for the instanton itself. The Yang-Mills action is invariant under the 15-parameter conformal and the 3-parameter 
global SU (2) groups. However, an instanton solution breaks 8 of these symmetries leading to 8 collective coordinates, 
roughly identified as 4 for its location, 1 for its size, and 3 for its isospin orientation, (see below for a more precise 
definition.) A monopole solution in each instanton background with these 8 coordinates held fixed further breaks 
some of the remaining 10 symmetries, leading to additional collective coordinates for the orientation of the loop. 

Next we will address the dependence of the loop solution on the loop radius R. Each solution is a stationary point 
of the gauge fixing functional under the constraint of fixed radius. To restrict further the MA projection, one may 
also impose the condition that G is a global minimum. We therefore need to study G(R), < R < oo, and determine 
if there is a true minimum for a finite non-vanishing value of loop radius R. We shall show that, for a large instanton 
with size p, (or equivalently small monopole loops with radii R << p), there is a very weak dependence of the gauge 
fixing functional G in R, so the scale of the loops is "nearly" a collective coordinate. This near zero mode leads to 
"marginal instability" for the formation of small monopole loops. 

A. 0(4) Invariance and Orientation of the Monopole loop 

The fact that we constructed the monopole loop solution lying in the 3-4 plane is purely for mathematical 
convenience; a larger family of loop solutions can be obtained by applying 0(4) Lorentz transformations. Since 
50(4) = SUl(2) x SUr(2), it is possible to define two sets of mutually commuting angular momentum generators 
by L a L = —{i/2)fi"x tl d v and L a R = — (i/2)rj^x^d u . (In terms of the conventional "rotation" and "boost" operators, 

J and K, one has Ll = [J — K]/2 and Lr = [J + K]/2.) Using the spinor basis, it is trivial to show that an 
instanton configuration is invariant under an arbitrary rotation in SUr(2), generated by Lr. An instanton is also 
invariant under an SUl(2) rotation, (generated by Ll), provided that a corresponding isospin rotation is performed 
simultaneously. 

If we consider the kinematics of the loop in the 3-4 plane, ignoring the instanton background for now, clearly 
rotations in the 1-2 and 3-4 plane leave it invariant. In fact it can be easily shown that these are the only invariances 
and that the remaining four-parameter coset space SUl(2)/Ul(1) x SUr(2)/Ur(1) rotates the loop to an arbitrary 
plane. Thus the loop breaks both left and right chiral SU(2) factors in the 0(4) group, analogous to the way the 
't Hooft-Polyakov monopole breaks the SU{2) isospin group. Since a MA gauge also fixes a direction in the isospin 
space by identifying T3, the only remaining symmetry transformations belong to SUr(2). Thus given a monopole 
loop solution, all other inequivalent solutions can be obtained by performing 0(4) rotations belonging to the quotient 
space SU R (2)/U R {1). 

We now provide a few details on these transformations. Any rotation UeSUii{2), expressed in the conformal 
coordinates of Eq.(|l3|), acts on x a ^ to give x 1 ^ = x a 0,U^,^. Let us see what happens to a circular loop of radius R 

centered in the 3-4 plane: u (cr) = xo{cr) + iyo{o) = and v (cr) = z a (a) + ito(a) = Re la / 2 , where < a < 4tt. First 
consider a rotation, i?3(A) = 

e i\L R = e ar 3 /2 ) m 

the U(l) subgroup of SU(2)r. It simultaneously rotates the 1-2 and 
3-4 planes by the same angle A/2, i.e., uo — * u(er) = uo(<r)e _lA / 2 and vo — > v(cr) = vo(cr)e~ lA / 2 . This clearly leaves 
a monopole loop lying in the 3-4 plane invariant. Next consider rotations, i?2(A) = c iXLr = e aAr2 / 2 . Again using 
conformal coordinates, one finds that 

u(o-) = cos(A/2)u (cr) - sin(A/2)v5(er) = - sin(A/2)i?e"' i ' T/2 
v(ct) = cos(A/2)v (cr) + sin(A/2)u5(cr) = cos(A/2)i?e JCT/2 . 



8 



The resulting loop has a circular projections onto the 1-2 and 3-4 planes with radii i?sin(A/2) and i?cos(A/2) 
respectively. In particular, for A = tt, it rotates a loop in the 3-4 plane to one lying in the 1-2 plane, as 
promised. Finally, consider a general SUr(2) rotation. Since it can be parameterized using Euler's angles as 
U(\3, A 2 , A3) = -R3(A3)i?2(A2)i?3(A 3 ), one finds that the resulting loop is given by 

u(o-) = -i?sin(A 2 /2)e-' i(<T+A3 -^ )/2 , v(a) ^ Rcos(X 2 /2)e l{a - X3 - y ^ /2 . 

For A3 7^ 0, although the projections onto the 1-2 and 3-4 planes remain circular, it is no longer circular for projections 
onto other planes, e.g., the 1-4 plane. Clearly, all distinct loops can be characterized by two independent angles, 
parameterized by A3 and A 2 . 

Given a fixed isospin orientation for an instanton, the average over the loop orientation tensor can now be 
found. For instance, for the standard isospin orientation given by Eq. (|^), the monopole loop solution lying in the 3-4 
plane is characterized by an anti-symmetric tensor = (6^.36^^ — 5^, 3^4). Averaging over SUr(2), one finds that 

< >= |[(^/i,3^i/,4 - £/i,3<W) - (<W^".2 - <W<W)]- ( 29 ) 

In a random dilute instanton gas, this implies a correlation between the monopole loops and the isospin orientation of 
the associated instanton. As discussed in the Conclusion this correlation may provide a signature of our mechanism 
for loop formation in the QCD vacuum. 



B. Limit of Small Monopole Loop 

Classically SU(2) gauge theory in four-dimension has no dimensionful parameter, so the instanton solution breaks 
scale invariance through the introduction of the width p. Consequently our monopole solutions depends on the 
dimensionlcss ratio p/R, and we may consider the solutions in the limit for small loops size, p/R — ► 00. Let us return 



to Eq. (27), where the width parameter p enters through the function f(x) = p 2 /[x 2 (x 2 + p 2 )]. After scaling both x 
and p by R, one finds that the equation is greatly simplified in the limit p/R — > 00, f(x) — » x~ 2 . 
In this limit, an exact monopole loop solution to Eq. (|27]) is 

0(x, 9) = (3 (x, 6) = 29- {9+ + 0_) + tt, (30) 

with 9± — tan~ 1 [u/(w ± R)], as defined for the earlier QED example with a monopole loop in the 3-4 plane. Note 
that /3o = on the u-axis, and it has a jump by tt on the v-axis at v — R, (/?o = —tt for < v < R and (3q = for 
R < v < 00). This solution is valid provided that < R << p. The Dirac sheet lies in the 3-4 plane bounded by the 
circle of radius R where the monopole current resides. 

In the limit of small monopole loop and for < x < < R, the function (3 ~ 29 — tt and therefore (3 is not "small" . 
This reflects the fact that O ~ (x) in the limit x — > 0, so it is singular at the origin. On the other hand, outside 
of the monopole loop radius, R < x < 00, (3 scales with R. Therefore in this "outer region", for R — > 0, (3 admits an 
expansion in (R 2 /x 2 ) as 

"»< 3 ~ - gjgagj) ~ " sm(2fl)( 5 ) + 0< £ ) ' (31) 

Consequently for fixed x and x =/= 0, the limit of small monopole loop is characterized by "small" gauge transformations. 
This allows us to examine the stability of monopole solutions by a linear analysis about (f> a (x) ~ 6 a ,3- 



C. Marginal Stability of the Small Monopole 



Our monopole solutions are stationary values of G constrained by the boundary condition to a fixed radius R. We 
therefore proceed to study the functional G(R) evaluated at the the loop solution in the range < R < 00. Near 
R — > 00, G(R) is monotonically increasing since G is divergent for an instanton in the non-singular gauge. For R 
small, a leading order monopole solution is known. Initially we had hoped that G(R) has a minimum for some fixed 
R > 0. However in spite of our best effort, variational calculations have so far led to results where G(R) is always 
monotonically increasing. This is also confirmed by very accurate numerical integration of our 2-d PDE's as described 
in Sec. 4. We are now convinced that the strong MA projection defined by the global minimum of the functional G 
is simply the singular gauge itself, which can be viewed as the limit of a monopole loop with radius shrinks to zero. 
Using our small radius solutions for (3, we are able to analyze the small R behavior for G(R). 
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We now consider in detail the region near to R = 0. Assuming that R = is a minimum, one would expect that 

G(R) ~ p 2 [So + Si(R/p) 2 + ], with d% > 0. It came initially as a surprise that we found Si = 0. This suggests 

the possibility of a zero mode in the stability equation around the singular gauge. 

Let us expand G[4>] to quadratic order in the neighborhood of <f)Q = (0, 0, 1), parameterized by </> ~ </>o + uJ, 

G[4>] ~ G[0 O ] + \ J d A xuo ■ Mlj, (32) 

where the constraint \<p\ 2 = 1 is realized to this order by having uj ■ cf>o — 0, i.e., w only has 2 transverse components. 
The stability of small oscillations is studied by finding eigenvalues of a hermitian operator, 



Muii = XiCOi, (33) 



where 



M = -d 2 - ~d x - - &f(x)T 3 L 3 L , (34) 

x x l 

with Ll as generators of SUl(2) and T of isospin rotations in the adjoint representation. 

Expanding ul in terms of these normalized eigenvectors, uj = ctiUJi, one has G[<f>] ~ G[</>o] + j Yli ^i\®i\ 2 - With 
R = as a minimum for G, stability requires all eigenvalues are non- negative, i.e., A^ > 0. One would also expect 
that an infinitesimal loop is "nucleated" along the direction uj of the eigenvector with the lowest eigenvalue, Ao, and 
ao = 0(R). It follows that Si oc Ao- The fact that our numerical treatment indicates that Si ~ as depicted in Fig. 
|3|, (described in Sec. 4) implies that Ao = 0, i.e., to quadratic order there is a zero mode. We now try to construct 
this zero mode. 

With uj in the 1-2 plane, T 3 takes on eigenvalues t% — ±1. Since M. commutes with L 2 L , L\, and T 3 , there are a 

family of eigenvalues {\{\ for each set of {1(1 + l),l^,t^}, I = 0,1,2, , — / < I3 < I, and they can be found by 

solving an ordinary differential equation, 

{-d x 2 - -D x + Ifcil _ 8t 3 hf(x)}V(x) = \V(x), 
X x z 

where ^ is the "radial" part of the eigenfunction. With f(x) > 0, the attractive case corresponds to ^3 = 1 and 
the lowest energy level is I = 1. We find that all eigenvalues are positive, except possibly one. This zero eigenvalue 
solution can be constructed by changing variable to z — 1/x 2 , and solving 

where V e ff(z) = 2/z 2 — 2/z 3 /(l/y / z) = 2p~ 2 /(z + p~ 2 ). The zero eigenvalue solution is 

*(z) = (l + ^)-(l + y)log(l + ^). (36) 

The x — > 00 limit is \&(x) — > as it must. But the x — > limit is ^(x) x~ 2 , which leads to logarithmic divergence 
in the norm at short distances. This divergence seems to render this "near" zero mode questionable. However, it 
should be recognized that our "small loop" solution (|3l]) is valid only for R < x < 00. A proper treatment of our 
small oscillations problem requires cutting out a small ball around the origin, with a radius tq = O(R). Moreover, 
because of the zero eigenvalue the divergence in the norm does not imply a divergence in G. We should in principle 
be able to solve this problem by find the lowest eigenvalue and then letting r go to zero. With respect to this cutoff, 
our "near" zero mode is physically meaningful. An alternative approach, which is under investigation, may be to use 
the more general MA gauge with a Higgs potential. This has a natural length scale l/m#, which probably provides 
the necessary cutoff at short distances. 

As one extends beyond the quadratic approximation, one can see that this zero mode leads to our "exact solution" 
in the special case of / = 1/x 2 . That is, the formal p — > 00 limit solution. In this sense, it is easy to see that the 
divergence should be cut off by p. For the exact solution, the eigenvalue is of the order (R/p) 2 log(p/R), which is 
consistent with our numerical calculation as exhibited in Fig. [| 
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FIG. 3. For A = 0.01, 0.10, oo, the change in the Gauge fixing functional G as a function of the monopole loop radius R. 



IV. NUMERICAL SOLUTIONS TO MA PROJECTION 



For a general value of the monopole loop radius R, we were unable to find an analytic solution. Consequently 
we have discrctizcd our PDE's and found numerical solutions. The need for a numerical integration method is not 
surprising since even the 't Hooft-Polyakov monopole has no known solution in closed form, except in the BPS limit. 

With the simplifying ansatz (p(f), we were able to reduce the problem from a 4-d to a 2-d set of PDE's allowing 
us to construct very accurate solutions on a 2-d grid. To check the validity of this ansatz, we also consider the 
standard 4-d hypercubic grid conventionally used in Monte Carlo studies of non-Abelian gauge field theory. An 
important advantage of the 4-d grid is the ability of making a global search for stationary points. However this 
numerical integration method should not be confused with Monte Carlo simulations. Here the grid is used merely to 
solve numerically the classical PDE's of Yang-Mills theory. As always for discretization methods, it is important to 
consider carefully errors arising from the grid spacing a and the volume of the box. Our analysis of these errors will 
shed some light on earlier investigations |15|-|l7j] of the Abelian projection in "cooled" instanton configurations. 



A. Single Instanton Case 

Once we have made our ansatz (p6[), /3(x,y, z,t) only depends on two "radial" co-ordinates: u = \j x 2 + y 2 and 
v = \J z 2 + t 2 . This feature also generalizes to the case of the MA projections (||) with a Higgs potential which allows 
the magnitude 4>(x,y,z,t) = |0| to fluctuate. Here our ansatz ( |26| ) implies that both P{xp) and 4>(Xfj_) depend only 
on u and v. The functional G takes the simple from 

G = 47r 2 J dudv uv [\M((3)<j? + i(<9 M 0) 2 + \\{<t? - l) 2 ] , (37) 

where cj){u,v) = \<f>\, {d^cj)) 2 — (d(f>/du) 2 + (d(j)/dv) 2 and 

M{8) = 8f(x 2 ) + (df3/du) 2 + (dp/dv) 2 + sin 2 /?(^ + 4) " 8/(x 2 )[sin 2 (3 - ± sin/? cos /3(- - -)]. 

U V ' U V 

The conventional MA gauge is given in the limit A — > oo. 

We introduce grids in u and v, mapping the infinite u — v plane to a finite region, for a range of values of R/ p and A 
(see Fig. ||). For all A (including the BPS limit A = 0), we found that the monopole loop solution exists and that the 
functional G increases as R A within errors. For all values of A the global minima appeared to be at the point where 
R — > and the instanton returned to the singular gauge. If we perturbed the background field slightly away from a 
pure instanton by changing the functional form, f(x) — p 2 j\x 2 (x 2 + p 2 )], one can easily find field configurations in the 
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same topological sector which stabilized the loop at a fix radius. This strongly suggests that quantum fluctuations or 
instanton interactions can cause monopole loop formation. 

To make a global search for minimal solutions for the MA projection and to verify our theoretical ansatz, we have 
resorted to a 4-d grid. A natural choice is the standard Wilson approach with link variables, 

U ll (x)=e ia M x ), (38) 

where a is the lattice spacing. On this grid the MA functional which leaves the Maximal Abelian subgroup (7(1) 
unconstrained is 

G = \ fl2 E^ 1 - \Tr$(x)U»{x)4>(x + li)U^ X y)}. (39) 

U,X 

In order to numerically approximate its value, we restrict the summation over a finite volume V around the origin 
with open boundary conditions on the surface. This is a good approximation for an isolated instanton, since the 
contributions to the sum for large x drops as 1 jx 1 . Furthermore to obtain a finite value for G in the infinite volume 
limit, we know that the gauge rotation must become a constant at infinity. The minimization of the restricted 
functional was done using the standard over-relaxation algorithm . 

The instantons are placed on the grid in the non-singular gauge and then rotated to singular form. Each link is 
approximated by its integral using the trapezoid rule. The resulting configurations gave topological charge and action 
deviating at worst by 10% from their continuum values, 1 and 2tt 2 respectively. 

We have investigated a range of different sizes of lattices and instantons. As expected |l5|Jl6fl , we have found 
monopole loop formation in the Maximal Abelian gauge. But our monopole radii do not scale with the instanton size. 
This is a clear signal that the stabilization of the monopole loop is a lattice artifact. Despite that, we find that once a 
loop is formed the gauge rotation accurately satisfies our theoretical ansatz ( pr^ ) that a — ip — ip and (3 depends only 
on u and v. 

A finite volume and finite spacing analysis has been done in order to further establish the fact that the monopole 
loops is a lattice artifact. We measure the change, 

AO/a, L) = , (40) 

instanton 

of the MA functional after gauge fixing and we demonstrate that it vanishes in the combined limit of infinite volume 
and zero lattice spacing. In order to show this, we have found a series of solutions for lattice sizes 20 4 , 22 4 , 24 4 , 26 4 , 
28 4 , 30 4 and 32 4 with a fixed instanton size. The infinite volume limit is obtained by doing a linear extrapolation 
in 1/L , where L = I/ 1 / 4 is the linear dimension of our lattice. Fig. [|(a) shows this extrapolation for p/a = 4. 
The second extrapolation to zero lattice spacing is done from these extrapolated values of A for various values of 
the instanton radius p. At infinite volume, the only scale in the problem is p, and it is therefore defines the lattice 
spacing, a ~ 1/p. Thus one takes the a — > limit by taking the p/a — ► oo limit. This extrapolation is done using 
instanton sizes p/a — 3, 4, 5, 6, 7. 

The results are shown in Fig. |](b). The extrapolated value of A(p/a,L) in the limit of a — > and L — > oo is 
clearly zero within the numerical error. Another indication that the global minimum of the MA functional is the 
instanton in the singular gauge is the fact that A is so small O(10 -3 ). This is also supported by our experience with 
the minimizing routine. Starting with the gauge background for an instanton in the singular gauge, it took very few 
iterations to converge which shows that the displacement of the minimum due to lattice artifacts is actually very 
small. 
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FIG. 4. The finite lattice spacing and finite volume dependence for the monopole loop solution on the lattice. (a) Typical 
Infinite volume extrapolation, (b) a — + extrapolation. 

B. Solutions for Instanton Pairs. 

We have also used our 4-d grid to begin to investigate solutions to the MA projection in the background of interacting 
instanton pairs. For example, we have begun to study instantons (I-I) pairs and instanton anti-instantons pairs (I-A) 
with various sizes, relative separations and relative isospin orientations. Already several general conclusions can be 
drawn. The I-I pairs appear to be very similar to the case of single instantons. Indeed for the 't Hooft ansatz, one 
can again show that the two instanton solutions (or indeed the multi-instanton solution), 

\-r = r a fj^djn(l + p\/x\ + p\lx\\ (41) 

already satisfies the MA gauge condition. On the 4-d grid our preliminary study indicates that the MA projection 
of the two instantons have small loops centered at each instanton which again appear to be due entirely to lattice 
artifacts. It is natural to ask whether there is a general property of all exact classical solutions that the global 
minimum of the MA projection has no monopole trajectories. 

On the other hand, we have computed the MA projection for I-A pair with isospins oriented in the most attractive 
configuration. In this case, contrary to the single instanton and I-I pair large loops are formed in the MA projection. 
At large separation, each instanton in the I-A pair has its own monopole loop as illustrated in Fig ^, but as the 
separation d is decreased the individual loops fuse at d/p ~ 1.88 to form a large loop which surrounds the I-A pair. 
These loops clearly scale with the size of the system and the reduction of the MA functional is clearly larger. Further 
more this reduction in G becomes larger as you approach the continuum limit. Although we have not yet finished 
the same detailed finite size and finite lattice spacing analysis as in the single instanton case, we have considerable 
evidence that this effect will survive in the continuum. 

In the dipole interactions between I-A pairs, we believe we are seeing the first semi-classical mechanism for nucleating 
monopole loops. As the instantons become denser, the monopole loops begin to percolate between the individual 
instantons. It is easy to imagine that in an instanton liquid there is a critical density at which the percolation clusters 
are infinite and the monopoles can be said to condense. This effect is an intriguing possibility for a semi-classical 
mechanism for confinement. Additional studies of monopole trajectories for interacting instantons and further analysis 
to support (or refute) this scenario are postponed to a future publication [i8fl . 
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FIG. 5. Monopole loops for I-A molecule with the most attractive orientation: (a) Loops stabilized at separation d/p = 1.: 
(b) Loops fusion at separation d/p = 1.38. 



V. CONCLUSIONS 



The main goal of this paper was to find the earliest point in the semi-classical instanton vacuum in which the 't 
Hooft Mandclstam monopole appears. We believe that this is the formation of small current loops centered at each 
instanton and anti-instanton. In the extreme limit of infinite separation, i.e., an isolated instanton, the loops shrink to 
zero. However to 0(R A logi?) in the radius there is a flat direction for loop formation in the gauge fixing functional. 
We have begun to investigate the effect of instanton anti-instanton interactions. Here it appears that when the relative 
isospin orientation of an instanton anti-instanton (I-A) pair is most attractive, the monopole loops are stabilized at a 
finite radius and as the pair moves closer together the individual loops fuse. 

There is much more to learn about the mechanism for loop stabilization and the effects of instanton interactions. 
For the most part we postpone this to a future publication [ fLs] , However it is worth posing some of the questions 
and extending several arguments touched on early. 

First we have noticed that the single instanton in the singular gauge already satisfies the differential form of the 
MA gauge and gives rise to a finite contribution to the gauge fixing functional G. In fact it is trivial to see that 
all multi-instanton configurations that satisfy the 't Hooft ansatz A® = Vj^dylniX + F(x 2 )) likewise satisfy the MA 
gauge with a finite contribution for G. Thus it is tempting to conjecture that any exact self-dual classical solution 
minimizes this functional without explicit monopole currents. We are investigating this conjecture further. In this 
scenario the essential mechanism for monopole loop formation would be the interaction terms between self-dual and 
anti-self-dual regions that act as the "domain walls" between I-A pairs. 

Next it is worth expanding a little on the relationship between the topological charge and the monopole charge for 
an isolated instanton. In the singular gauge, when we write the topological charge, Q = (1/167T 2 ) / d 4 ccTr[F Ali ,-F Ml/ ], in 
terms of the gauge variant current, Q = J d xd^K^, one must exclude singular regions. Indeed via Gauss' theorem, 
the flux into these excluded singularities gives the net topological charge. We have studied carefully how this theorem 
is satisfied in the presence of our monopole loop. In the singular gauge, we may write the instanton as 



A„ 



i x 2 + p 2 



(42) 



As we shrink a small sphere toward the origin x 2 
singular part of the current, 

1 



the surface area is 0(8 3 ). However one sees that the 



K n = g - 2 t/ivpxTr \A V (dpA\ - i^A p A\)], 



(43) 



diverges as 0(1/S 3 ) and that the contribution to Q comes entirely from the pure gauge piece, M M . On the other hand 
as we rotate the field configuration into a monopole loop, 
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(44) 



the singularity is spread out to the loop at radius x 2 — R 2 . Now the toroidal surface area around the loop is 
0(RS 2 ), but again a careful analysis shows that the only divergent piece of 0(l/5 2 ) comes from M^. The argument 
is completely general for a monopole loop gauge and it requires the magnetic charge to be quantized to match the 
topological charge. This shows that the monopole loop has a very natural fit to the instanton. 

Finally we have found that when an isospin orientation of an isolated instanton is fixed, the monopole loops are 
constrained to a portion of the full 0(4) group. Again there is a rather remarkable "coincidence" this time between 
the symmetries of the monopole loop and the instanton. A single loop "kinematically" breaks the Lorentz group 
SUl(2) x SUr(2) into a factor of coset spaces SUl(2)/Ul(1) x SUr(2)/Ur(1) in the same sense that the 't Hooft- 
Polyakov monopole brakes the SU(2) isospin symmetry down to SU(2)/U(1). Then if we fix the isospin axis of the 
instanton (anti- instanton) and choose the MA gauge relative to T3, the loop is only re-oriented by the right (or left) 
coset. This implies a correlation between the plane of the loop and the isospin orientation, which can be tested in 
typical background configurations of instantons generated in a Monte Carlo simulation. In this way we can determine 
whether or not our conjecture that these monopole loops are important configurations for the full quantum theory is 
correct. 

In conclusion, although there are many more details worth considering there is a remarkable coincidence between 
the form of an instanton and its monopole loop in the MA projection. This is reflected in topological, symmetry and 
stability terms. This leads us to see the instanton in a new light as the "seed" for the formation of monopole loops. 
The dynamical implications are much more difficult, but it appears that I- A interactions may play a crucial role and 
the large "entropy" of monopole loops percolating between near by instantons suggest a promising direction for future 
research on electric confinement. 
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APPENDIX A: GAUGE COVARIANT FORMULATION OF THE ABELIAN PROJECTION 

The usual way to study the MA gauge is to applying the gauge transformation to the field, 

AJx) -» A^(x) = (l(x)AJx)nUx) + -n(x)dMUx) . (Al) 

and find that rotation f2(x) which satisfies a gauge condition. Alternatively, one may introduce an auxiliary adjoint 
Higgs-like field, $ = <f> ■ t , expressing the entire problem in a gauge invariant form. We refer to this two methods as 
"active" and "passive" respectively. In this appendix, we collect together basic formalism for each form. 

In the passive description, the MA gauge is cast in a gauge covariant form familiar to the Higgs |^5j model, with 
the entire formalism re-expressed exactly in the form used to construct the 't Hooft-Polyakov monopole. This 
covariant formulation both simplifies our analysis and emphasizes the important physical point that the Maximal 
Abelian projection needs not be viewed as a gauge fixing prescription but rather as a way to define magnetic degrees 
of freedom. We prefer the latter interpretation!] 



3 There is also a considerable literature |2^,^], which goes on to try to construct an appropriate monopole condensate order 
parameter in the (7(1) sector. Our field <£> is not this object, since it, like the Higgs field, lives in the coset space. 
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1. Active View of MA Gauge and Magnetic Current 



In the active form, the resulting Abelian field a fl (x) = Tr[rsA^(x)], (or more properly its derivatives), which 
depends on the choice on the non-Abelian background field, can have singularities leading to a non-zero magnetic 
current. Introducing the notation, A^x) = fl(x)A ll (x)il ji (x) and M M (x) = (l/ie)fl(x)d^' f (x), the transformed field 
becomes 

A%(x) = A^x) + M^x). (A2) 

This splits the Abelian field into two components, a fl (x) — A^(x) + M^(x). The first term, A^ 7 for our problem, 
comes from a direct rotation of the instanton field. The second "induced" term, M^, can contain monopoles as its 
source when appropriate conditions are met as we demonstrate next. 
The Abelian field strength, = d^a^ — d„a^, is given by 

/„„ = (n^nt) 3 - ie[n{A^ + ^)n\n{A v + ^d u )n% (A3) 

where (• • -) 3 = Tr[r 3 ■ ■ ■]. It again is split into two pieces, 

/„„ - (d^Al - d v AD + {d„Ml - d v Ml). (A4) 

Since the dual of the first combination is obviously divcrgcntlcss, only the second term contributes to a non- vanishing 
magnetic current, 

fc M (z) = ^e^ pa Tr(T 3 d4M p (x),M a (x)]), (A5) 
where = j^d^f u ^, and f ufl = T^tpupo fpa ■ 



2. Gauge Invariant View of MA Projection 

We now reformulate the MA projection in passive form. Let us begin by noting that the functional, 

G=\ J {Al + iAl)(Al-iAl)d 4 x (A6) 

is just the mass term in the broken phase of an SU (2) Georgi-Glashow model. This suggests a change of variables 
from Cl(x) to 

$(jc) - n\x) T3 n(x). (A7) 
The functional G now takes form of action for a Higgs field, 

G = i J <fx{\[D„{A)$Y + V(?)}, (A8) 

where $ = D^4> a = + ee abc A^cff , and the potential, 

V(P) = *(& 1), (A9) 

involves a Lagrange multiplier a in order to maintain the constraint <fi • <fi — 1 ■ We also suggest a generalization of the 
MA projection to include the standard quartic Higgs potential, 

V^^^-v 2 ) 2 , (A10) 

with $(x) now given by (f>(x) ■ t/\4>(x)\. This more general from we will refer to as the "Higgs MA projection". 
Thus the MA gauge is precisely the same as minimizing the action for a non-dynamical (or auxiliary) Higgs-likc field 
in a background Yang-Mills theory. The limit A — > oo at fixed VEV v results in the usual MA projection. This 
is the limit which takes the Higgs mass nig to infinity. This formulation corresponds to a passive description of a 
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gauge transformation. In the active description, rotation Q is applied directly to the gauge field in Eq. (||) so that 
(f> — ► (0, 0, 1). Clearly, the passive and the active descriptions are equivalent: each specifics a gauge transformation Q, 
up to the U(l) subgroup. Ho weve r the passive description is more "natural" since 4> lives in the coset space. 
The Abelian field strength (A3) now takes a manifestly gauge invariant form, 



^abcta 



D^°D v <j) c , 



Uu = - iee a 

and the magnetic current is the well known conserved topological current, 

1 



-e^ pr7 e abc d^ a d P (f> b d^ 
8ne 



(All) 



(A12) 



The essential identity in deriving these expressions from our earlier equations is 9^$ = i[f^A/ M f2, $]. An alternative 
form for f^v follows directly from Eq. ( ]A4| ) , 

Uu = d^Al - d„Al - iee abc rd^ b dj c . (A13) 

This approach to MA projection not only provides a gauge invariant treatment, but also allows a topological interpre- 
tation for the possible presence of magnetic sources. For example, as in the case of the 't Hooft-Polyakov monopole, 
the presence of a monopole charge can be understood to be due to a non-trivial homotopy U2(SU(2)/U(1)) = Z. 



3. Higgs MA Projection for Instanton 

For general reference to our analysis, we final give the full set of differential equation for the Higgs MA projection, 

D,j,{A) 2 4>- W^V{4>) = . (A14) 
The two independent PDE's for the tangential components, now take the form 



(j> d 2 a + 2d^d p a + 2(f) cot/3 (d^a)(d^/3) 



5 2 /3 + 2d^(3 - ±0 sin2/3 (d^a) 2 = 2[0sin/3(0 • A p ) {d„a) - (a ■ 1^(0^)], 



(A15) 
(A16) 



by projecting Eq. (11) onto the a{x) and P(x) axes respectively. In addition there is a new equation for the radial 
mode, 



dlcf> = -{2[ S in/3(/3 • A^a) + (a ■ A^(3] + - U 2 ]}^ 



(A17) 



Under our ansatz, a — ip — ip, with (3(u,v) and <p(u,v) only functions of u and v, Eq. ( A1E ) is still satisfied 
automatically. Now Eq. (A.16) becomes 



%P + 2(d u <t>d u f3 + d v cj ) d v p)-{ 



^ir)0sin(2/?) 

lU V 



V{x){4>[(- 



-) sin 2 (3 — sin 2/3] + [vd u (f> — ud v cj)]}, 



and the the radial equation (A17) becomes 



(A18) 



9> = -{/(x)[sin2/3(- 



-) + 4sin^ (3 + 2{ud v [3 - vd u f3)] + \V\F) - U 2 ]}<p 



(A19) 



Both analytical and numerical properties of these equations are discussed in the text. There appears to be remarkably 
little dependence of our monopole loop solution on A from A = (the BPS limit) to A = oo (the standard MA 
projection). 
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4. Abelian Projected Theory 



In the MA gauge, one is left with an intermediate description of a U(l) gauge theory, 

C{a»,A±) = \fi v + (d^A+^A^ + eHA+A- - A+A~)A+A~. (A20) 

interacting exclusively through charged vectors A* = A 1 ^ ± iA't. The U(l) covariant derivative is d^ = 9 M + iea^. 
In addition to the standard Wilson loop W(C) — TrVc^ $ dx^A^ ^ Qne ma y m ^ roc j uce fae Abelian Wilson loop, 

W ab {C) = e */ dx M a M (A21) 
which in turn can be split up into "photon" and "monopole" factors 

W ab {C) =exp[ij dx^Al] exp[ie J J da^e abc ^ a d^ b d,4> c ]. (A22) 

The second integral is taken over any surface whose boundary is given by the loop C. The near "saturation" of the 
non-Abelian string tension by the Abelian part in the MA gauge, <r a b ~ .92a, and the Abelian part by the monopole 
contribution, o monopo i e ~ .95cr a h, observed in lattice simulations pgj is referred to as Abelian Dominance. If this 
survives the continuum limit this may provide the dynamical link between monopole configuration and confinement. 

To understand the role that Abelian dominance might play in the continuum theory, it is useful to look at the very 
interesting form of the Wilson loop suggested by Diakonov and Petrov [^9) , 

W(C) = J V$(x) exp[i§ J dxj ■ A M ] exp[ J da^e abc 4> a dj b dj c }. (A23) 

A comparison with the Abelian Wilson loop in the MA gauge shows that this has exactly the same form with the 
additional step that one must average over all gauge transformation in the SU(2)/U(1) coset. Consequently Abelian 
dominance is the statement that in the true quantum vacuum, the contributions to the <f> average is approximated 
by the MA projection. Thus one strategy to proving confinement is to first establish Abelian dominance (or more 
precisely the inequality W(C) < W ab (C) for large loops) and then to demonstrate that monopole condensation occurs 
forcing an area law for the Abelian loop. 
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